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We derive a new perturbation scheme for treating the large d limit of lattice models at arbitrary 
filling. The results are compared with exact diagonalization data for the Hubbard model and found 
to be in good agreement. 
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I. INTRODUCTION 

In recent years there has been a renewed interest in the 
study of doped transition metal oxides like Lai^xTixO^. 
These materials exhibit interesting phenenoma like the 
correlation induced metal insulator transition. Although 
there are several experimental data available right now 
|^|-^ it is still quite difficult to tackle these substances 
theoretically. Realistic models have to take into account 
several bands and are to be explored at finite doping. 

The most promising way towards a theoretical descrip- 
tion, perhaps, is the limit of large spatial dimensions 
which defines a dynamical mean field theory for the prob- 
lem. This limit can be mapped onto an impurity model 
together with a selfconsistency condition which is char- 
acteristic for the specific model under consideration 
The mapping allows to apply several numerical and ana- 
lytical techniques which have been developed to analyse 
impurity models over the years. There are different ap- 
proaches which have been used for this purpose: qual- 
itative analysis of the mean field equations j5|, quan- 
tum Montecarlo methods 1^-^, iterative perturbation 
theory (IPT) |^,^, exact diagonalization methods 
and the projective self-consistent method, a renormaliza- 
tion techique . However, each of these methods has 
its shortfalls. While quantum Montecarlo calculations 
are not applicable in the zero temperature limit, the ex- 
act diagonalization methods and the the projective self 
consistent method yield only a discrete number of pols 
for the density of states. Moreover, the computational 
requirements of the exact diagonalization and the quan- 
tum Montecarlo methods are such that they can only 
be implemented for the simplest hamiltonians. To carry 
out realistic calculations it is necessary to have an accu- 
rate but fast algorithm for solving the Impurity model. 
In this context iterative perturbation theory has turned 
out to be a useful and reliable tool for the case of half 
filling jl^JT^. However, for finite doping the naive ex- 
tension of the IPT scheme is known to give unphysical 
results. There is still no method which can be applied 
away from half filling and which at the same time is pow- 
erful enough to treat more complicated models envolving 
several bands. 

The aim of this paper is to close this gap by introducing 
a new iterative perturbation scheme which is applicable 
at arbitrary filling. For simplicity, we treat the single 
band case here. But we believe that the ideas can be 
generalized to more complicated models involving several 
bands. 

The (asymmetric) Anderson impurity model 
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describes an impurity {fa) coupled to a bath of conduc- 
tion electrons (cfeo-)- The hybridization function is given 

by lS.{uj) :— J2^. ^^^^ . Once a solution is known for arbi- 
trary parameters a large number of lattice models can be 
solved by iteration. An example is the Hubbard hamil- 
tonian: 
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which can actually serve as an effective hamiltonian for 
the description of doped transition metal oxides . On 
a Bethe lattice with infinite coordination number z the 
Hubbard model is connected to the impurity model by 
the following selfconsistency condition: 



(3) 



and 6/ = — /i. The mapping requires that the propagator 
of the lattice problem is given by the impurity Green 
function (G — Gf). Below we set D := 2t — 1. 

In the next section we will derive the perturbation 
scheme for the impurity model. Afterwards the scheme is 



applied to the Hubbard model (section [II ) . Some results 



for the doped system are presented and the accuracy of 
our scheme is discussed. We conclude with a summary 
and an outlook on further extensions (section 



II. DERIVATION OF THE APPROXIMATION 
SCHEME FOR THE SINGLE IMPURITY MODEL 

In this section, we derive the approximation scheme 
which, given the hybridization function A(lu) and the 
impurity level e/, provides a solution of model (|l|). For 
simplicity, we assume that there is no magnetic symme- 
try breaking (n^ — — n). We also restrict us to zero 
temperature. The procedure is an extension of the ordi- 
nary IPT scheme to finite doping. The success of IPT 
at half filling can be explained by that it becomes exact 
not only in the weak but also in the strong coupling limit 
|]l3| . Moreover, this approach captures the right low and 
high frequency behavior so that we are dealing with an 
interpolation scheme between correct limits. 

The idea of our approach is to construct a self energy 
expression which retains these features at arbitrary dop- 
ing and reduces at half filling to the ordinary IPT result. 

Ordinary IPT approximates the self energy by its sec- 
ond order contribution: S(cij) w Uu + 'Eq (lu) where 
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with — ilmGo. Here, the (advanced) Green function 
Go(a;) is defined by 
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Go(^) := 



1 

uj + fLo — A{uj) 



(5) 



The parameter /io is given by — e/ — Un. In particular it 
vanishes at half filling. The full Green function follows 
from 
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Go ^ - Ao - £/ - S(t^) 



(6) 



To ensure the correctness of this approximation scheme 
in different limits, we modify the self energy functional 
as well as the definition of the parameter fiQ. 

We start with an ansatz for the self energy: 



(7) 
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Here, Sp (^) is the normal second order contribution 
defined in equation (^). We determine the parameter A 
from the condition that the self energy has the exact be- 
havior at high frequencies. Afterwards, B is determined 
from the atomic limit. 

The leading behavior for large w can be obtained by 
expanding the Green function into a continuous fraction 

— L). Here, ai marks 



H: G/(fc,c.) = l/(c.-e/-ai~2^ 
the ith order moment of the density of states. One can 
compute these quantities by evaluating a commutator 
(see [0). We obtain for our model G/(fc = kfermi,^^) — 
l/(uj — ej — Un— ^ ) ■ The leading term of the self 

energy is therefore given by 



^Un + U^n{l-n) - + 




(8) 



Here, n is the physical particle number given by n = 
duj ImG f (uj) . has to be compared with the large 
frequency limit of (M): 



(9) 



where uq is a fictitious particle number determined from 
Go (i. e. no = J°^dulmGoH). From (0), (|), and (|) 
we conclude 
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Chosing A in this way guaranties that our self energy 
is correct to order —. It should be noted from the con- 
tinuous fraction considered above that consequently the 
moments of the density of states up to second order are 
reproduced exactly. 

Next, we have to fix B. The exact impurity Green 
function for Vk ^ is given by ll^] 
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This can be written as Gf{uj) — 1/{lu — ef — T^atomici^)) 
where 



Satormc(w) =^ Un + ■ 



n(l - n) C/2 



(12) 
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This expression is to be compared with the atomic limit 



of our ansatz (^. Since Tr^' (oj) 
B = 



U rto(l-Tio) 
LO+flo — iri ' 
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we obtain 
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no(l - no) C/2 
Thus, the final result for our interpolating self energy 
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Yet, jlo is still a free parameter. We fix it imposing the 
Friedel sum rule fB: 
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This statement, which is equivalent to the Luttingcr the- 

orem ^ tlGf{^f-^i^ = 0, should be viewed 

as a condition on the zero frequency value of the self en- 
ergy to obtain the correct low energy behavior. The use 
of the Friedel sum rule is the main difference to an ear- 
lier approximation scheme and is essential to obtain 
a good agreement with the exact diagonalization method. 

So far, we considered three different limits: strong cou- 
pling, zero frequency and large frequency. It remains to 
check the weak coupling limit. Taking into account that 
n = no and jlo = — £/ for C/ = 0, it follows that (Q) is 
indeed exact to order U^ . 

The actual solution of the impurity model is deter- 
mined by a pair [jlo, n) which satisfies equations (Q), 
(^), (^), (^4|), and (|l5|). For the numerical implementa- 
tion Broyden's method a generalization of Newton's 
method, has turned out to be very powerful. Defining 
two functions /i(n, jlo) := n—J^^ dci;ImG/[n, /io](i^) and 
/2(n, jlo) := npriedeiVi, jlo] the impurity problem can 
be solved by searching for the zeros of /i and /2 (n Friedel 
is the particle number determined from the Friedel sum 
rule). The algorithm is very efficient as in most cases a 
solution is found within 4 to 10 iterations. 



III. APPLICATION TO THE HUBBARD MODEL 

After treating the Anderson impurity model, we now 
apply the perturbation scheme to the solution of the Hub- 
bard model. Starting with a guess for /S.{uj) one can solve 
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the impurity model using the scheme described above. 
This yields a propagator Gf = G, which can be used to 
determine a new hybridization function A(a;) according 
to (^). The iteration is continued until convergence is 
attained. It is most accurate to perform the calculation 
first on the imaginary axis. Once the constants A and 
B in the interpolating self energy are determined in this 
way, they can be used to perform the iteration on the 
real axis. 

In the case of the Hubbard model, the Luttingcr theo- 
rem takes the simple form pl| 



Im G(C0) 



^0 = M - Si„i(a) = 0), /io:=/-i|(7= 



(16) 



This can be used to simplify the selfconsistency procedure 
if pLf) rather than /i is fixed. Starting with a guess for G 
and /i, one can compute Go, n, and tlq. Afterwards (|l^ ) 
yields Sj„t(a;) and a new ^ is obtained from (|l6|). 

To illustrate the accuracy of our method we compare 
it with results obtained using the exact diagonalization 
algorithm of Caffarel and Krauth . Both methods are 
in close agreement when used on the imaginary axis (see 
figure |l|) . The real advantage of our perturbation scheme 
compared to the exact diagonalization is disclosed when 
we display the spectral functions obtained by these two 
methods on the real axis (figure 0). 
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FIG. 1. lmG{iuo) at T = for 17 = 4D and hole doping 
5 = 0.14: iterative perturbation scheme (full line) vs. exact 
diagonalization (dashed line) 
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FIG. 2. lmG{Lu) at T = for [/ = 4D and hole doping 
5 = 0.14: iterative perturbation scheme (full line) vs. exact 
diagonalization (dashed line) 

It is clear that the exact diagonalization is doing its 
best in producing the correct spectral distribution. But it 
is unable to give a smooth density of states. Instead sev- 
eral sharp structures occur as a consequence of treating 
only a finite number of orbitals in the Anderson model. 

Figure || shows the evolution of the spectral density 
of the doped Mott insulator {U = AD) with increasing 
hole doping 5. The qualitative features are those ex- 
pected from the spectra of the single impurity j5| and are 
in agreement with the quantum Montecarlo calculations 
p2[ . For small doping there is a clear resonance peak at 
the fermi level. As 5 is increased, the peak broadens and 
is shifted through the lower Hubbard band. At the same 
time the weight of the upper band decreases. 

The most striking feature of the evolution of the spec- 
tral density as a function of doping is the finite shift of 
the Kondo resonance from the insulating band edge as 
the doping goes to zero. It was demonstrated analytically 
that this is a genuine property of the exact solution of 
the Hubbard model in infinite dimensions using the pro- 
jective self-consistent method and is one of the most 
striking properties of the Hubbard model in large dimen- 
sions. This feature did not appear in the earlier studies 
of Hubbard model in large dimensions using Montecarlo 
techniques |22j at higher temperatures, and is also not 
easily seen in exact diagonalization algorithms ll2]. 
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FIG. 3. Evolution of the spectral function for U = AD and 
T = with increasing hole doping S 



IV. DISCUSSION AND OUTLOOK 

In this paper we introduced a new perturbation scheme 
for the solution of lattice models away from half filling. 
The basic idea is to construct an expression for the self 
energy which interpolates between correct limits. In the 
weak coupling limit our approximate self energy is exact 
to order U"^, and it is also exact in the atomic limit. The 
proper low frequency behavior is ensured by the Friedel 
sum rule (or, equivalently, the Luttinger theorem). This 
is important to obtain the right low energy features in 
the spectral density. The overall distribution of the den- 
sity of states on the other hand is determined by the 
spectral moments, which arc reproduced exactly up to 
second order by satisfying the proper large frequency be- 
havior. In the light of these features it might not be too 
astonishing that we obtain a good agreement with the 
exact diagonalization method. 

Since the algorithm decribed here is accurate and very 
fast (a typical run to solve the Hubbard model takes 60 
seconds on a DEC alpha station 200 4/233) it has a wide 
range of applications. Two examples that come to mind 
are the effects of disorder on the Hubbard model away 
from half filling and the study of realistic models with 
orbital degeneracy. The latter is very important to make 



contact with realistic three dimensional transition metal 
oxides. 
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